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1 . Int r oduc t ion 


Fluid-dynamic forces arising from nonuniform pressure .patterns in labyrinth seal 
glands are known to be potentially destabilizing in high power turbomachinery. A 
well documented case in point is that of the space Shuttle Main Engine turbopumps 
(ref. 1), and other examples can be found in the literature, as for instance in the 
recent review of Ehrich (ref. 2) and Childs and Ehrich (ref. 3). Seal forces are al- 
so an important factor for the stability of shrouded turbines, acting in that case in 
conjunction with the effects of blade-tip clearance variations (refs. 4,5). 

The basic mechanisms which produce the uneven pressure distribution in a laby- 
rinth have been qualitatively or seml-qualitatively discussed in many references 
(refs. 6,7). In most instances, the basic agent is found to be flow swirl in the 
glands, either from pre-swirl (as in the case of turbine shrouds) or from frictional 
interaction with the rotating shaft (as in multichamber jet engine seals). Quantita- 
tive modeling of these forces has also been reported by a number of authors (refs. 8, 
9,10), using lumped-parameter models for each gland. These models yield in general 
predictions of the direct and cross-wise stiffnesses and damping coefficients for 
small shaft displacements, and are useful for linear stability analyses. Non-linear 
predictions for fully developed unstable operation are less advanced. 

Fairly extensive data also exist on the stiffness factors of seals of various geo- 
metries (refs. 4,5,11,12,13). These have been generally obtained in rotary rigs with 
adjustable shaft eccentricity. Much less satisfactory is the situation with respect 
to data on damping coefficients due to labyrinth seals, since these require dynamic 
measurements on either vibrating shafts, or shafts fitted with adjustable whirl mech- 
anisms. Yet these data are almost as essential as those on stiffnesses, since the 
corresponding induced forces are of the same order. Ref. (13) reports damping data 
for non-rotating shafts. 

In this paper we report on the development of a linear model for the prediction 
of labyrinth seal forces and on its comparison to available stiffness data. We also 
present a discussion of the relevance of fluid damping forces and report on the pre- 
liminary stages of a program to obtain data on these forces. 

2. Model Formulation 


The model is very similar in its main outline to those of Kostyuk (ref. 8) and 
Iwatsubo (ref. 9). It describes the flow of an ideal gas through the seal chambers, 
assuming largely constant temperature, but allowing for isentropic acceleration to- 
wards the narrow gaps and also for isentropic azimuthal flow redistribution in each 
chamber. Each chamber is assigned a pressure P. and azimuthal velocity c^, and these 
quantities are governed by equations of mass anS azimuthal momentum conservation, 
written in integral form. 

The axial flow rate q through each seal throttling is approximated by a commonly 
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used expression, which basically derives from Bernoulli’s equation with a density 
halfway between those of the two adjacent chambers. Per unit length in the azimuthal 
direction, this gives 
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where 6^ is the width of the narrow gap (see fig. 1) and ]i^ = c^ 3^ is the product of 
the usual contraction coefficient c times a carryover factor 3 i to account for 
nonzero upstream axial velocity and nonzero pressure recovery in the downstream cham- 
ber. Eq. (1) is assumed to apply locally at each time t and azimuth (f). This semi- 
incompressible approximation is known to be reasonable up to gap Mach numbers of 
about 0.5; however, the last gap or two of a labyrinth with a high overall pressure 
ratio may be above that M ach number, and, in particular, the last chamber may choke. 
We have partially accommodated this effect by retaining Eq. (1) throughout, but re- 
placing it by a choked-flow expression in the last chamber only if the first approxi- 
mation indicates sonic or supersonic conditions there. 

With reference to the geometry of fig. 1, the governing equations within each 
chamber are 
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These equations are first linearized about a condition of zero eccentricity. The 
zero’ th order approximation provides a basic flow rate q and pressure and azimuthal 
velocity distributions P* , c* , (Appendix 1) . The first approximation then provides 
linear equations for the perturbations, defined by 

Pi = p* (1 + = q* (1 + (i + (4) 

where rj^ are functions of t and w == R^(|) . The right hand sides of these 

equations are determined by an assumed eccentric motion of the shaft, whose center 
follows an elliptic path 

X = r cos flt ; y = r sin ^^t (5) 

C l c 2 

where H is the shaft vibration frequency, closely identified with one of its natural 
frequencies. 

The details of the analysis are given in references 14 and 15. For convenient 
solution, the perturbations (for a stationary oscillation) are expressed in the form 

5 - R [e “^(5 sin (}) + 5 cos 4>)1 » etc. (6) 

e s c 
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where R indicates the real part, and 5 » C itJ- general complex numbers. After 

elimination of , the following system of^ perturbation equations is obtained: 
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The expressions for the coefficients A. , B. , , Z. are given in Appendix 2. 

For a seal with K chambers, i would range from i = 1 to K. The calculations reported 
here have assumed uniform inlet and exit conditions. 

it 

The forces are then obtained by integration of the perturbation pressures 
around the periphery of each chamber, followed by summation for all chambers, Sxnce 
our attention will be restricted to the practical case of circular whirl (ri= Xz~ r) , 
it is advantageous to project the forces in the directions towards the instantaneous 
minimum gap (1*^) and 90° ahead of it in the whirl direction (F ) , as shown in fig. 2. 
We obtain ^ 

2tt ^ 

F- = -R £. / cos - fit) d(}) (8a) 
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X 0 

2ir ^ 

F = -R £. / P.5* sin (({> - Qt) d(j) (8b) 

qi s X Q X X 

These components can be expected to be time^invariant for a symmetric shaft in 
circular whirl. For this case, the form of the system of equations (7) indicates 
that we must have 
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and the system (7) reduces to 
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Also, the complex force vector + i F is found to be given simply by 

X ^i 


+ i F = 7T R P. (“ i ^ ) £. 

s X X 
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whereas, if the forces F and F along fixed directions Ox, Oy are desired, they 
are given by ^ ^ 


F 4* i F = (F^ + 1 F ) e 
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3. The stiffness and damping coefficients 

Most data reported to date on seal forces refer to situations with a static off- 
set, i-e,, with fl = 0. Within a linear approximation, it is then unambiguous to de- 
fine the direct and crosswise stiffness K K , as - F^/r and - F /r , respectively. 
In practice, however, unstable whirl is o^servM to occur at or neaS (one of) the 
shaft natural frequencies, usually at the first one, Q (and it first shows up at 
rotational frequencies O) of the order of twice this natural frequency) . For purposes 
of dynamic modeling, then, it is of interest to calculate or measure the forces for 
Q ~ , since the fluid disturbances are expected to translate into relatively small 

real and imaginary departures of Q from its basic resonance value. A common expres- 
sion for the disturbance forces in terms of the x , y deviations of the shaft center 
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where the stiffnesses K., and damping factors C, . are taken to be constant. This 
amounts to an approximation in which terms proportional to x , x , etc., are all 
neglected. It is easy to show that the corresponding approxxmatxon in terms of Q is 
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one where the actual functions F (fi) , F (ft) are replaced by their tangents at the fre- 
quency of interest: 


F , (ft) F , + F , (ft - ft ) 
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and then the K and C coefficients are given by 
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Thus, determination (either analytical or experimental) of F , F at two whirl- 
ing frequencies near Q is sufficient to extract the K and C coefricieSts in this 
formulation. 


An alternative formulation is often found in the literature (refs. 8,9,15) and 
is, in fact, the basis for the results presented here. Analytically, it consists of 
solving the system of equations (7) for V 2 = 0, i.e., for a linear vibratory shaft 
motion, and adscribing to^-K x , K x those forces F , F in time phase with x 
while attributing to -G x -C x^^ ^ those in quadrature ^ with x . It can be^ 
shown that this leads to the definitions 
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A geometrical interpretation of the difference between equations (15) and (16) 
is shown in fig. 3. Experimental or analytical determination of the set (16) of co- 
efficients requires data on F^ and F^ at both Q and -Q. 

The KVs and C's given by (15) and (16) coincide only if F^ and F are linear 
functions of Q. Since either definition may be used in the ^literature 

with little elaboration, it would be of interest to study the extent to which this 
leads to numerical differences; pending this, we will in this study adopt the defini- 
tions (16). An example for a single-chamber seal is shown in Appendix 3. From this 
limited evidence, it appears that the important coefficients K , G are about the 
same in both definitions. 
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4. Comparison to Literature Data 


4.1 Data Used 


Benckert and Wachter have published an extensive set of data (refs. 11,16) for 
multichamber labyrinth seals of simple "straight*' or "full'' types (fig. 4). The 
data were taken in a static-offset rig operating on pressurized air, and induced 
forces were obtained by integration of measured azimuthal pressure variations on a 
number of seal cavities. Labyrinths with up to 23 chambers were used. The experi- 
ments allowed variation of shaft speed, (O, overall pressure ratio P /P , rotor 
eccentricity r, seal geometry (6^ , , h. , Pig. 1) , number of ° chambers 

and entry swirl c . The seal flow rate q was measured and an averaged carry- 
over factor y waS deduced from these data and reported in a number of instances. In 
our calculations we used these "measured" factors when available directly; in other 
cases, we adopted values measured for chambers of the same geometry, or, for the 
"full" type of seal, where little carryover is expected we used y = 1. The contrac- 
tion coefficient c was taken as a function of Reynolds number and strip geometry as 
given by Vermes (ref. 17) (fig. 5). 

Data of Brown and Leong (ref. 18) were also used for validation of the undis- 
turbed flow predictions. 

4 . 2 Undis turbed Flow Parameters 


Figure 6 compares Brown and Leong 's data on the axial pressure distribution in 
an 11-chamber test seal with our calculated undisturbed pressure distribution. There 
is good agreement except for the sharp pressure drop shown by the data between the 
inlet and the first chamber. This is probably a reflection of a reduced carry-over 
factor on the first strip; the calculation used a constant y (the value used is ir- 
relevant to the comparison) . 

Benckert and Wachter reported for one particular case the axial variation of 

azimuthal velocities c* . This was for a 23-chamber seal with c* = 40 m/sec, 

P /P = .66, R = 0.15^, 0) - 1000 rad/sec, ri = 0.25 mm, Xz = 0, 6^ = 0.5 mm, 
a o s 

= 4 mm and h^ = 6 mm. The data are shown in figure 7, together with the code pre- 
dictions. The good agreement shown is important for the prediction of disturbance 
side forces, which depend critically on swirl velocities. These results appear to 
validate the formulation used for the friction factors between the fluid and the 
stator and rotor surfaces (turbulent pipe flow formulae with a standard correction 
for "pipe" curvature). 


4.3 Stiffness Coefficients Without Shaft Rotation 


for a number of cases from Benckert 


The cross-spring coefficients K = - K 
and Wachter 's tests with a non-rotating shaff^were calculated and the results are 
summarized in Table 2 and figure 8. The key in Table 2 describing the test para- 
meters is explained in Table 1. 


ri 


The eccentricity ri used in the tests was 0.15 mm, except for Run 17, which had 
= 0.25 mm. 
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Table 1 . Key for Table 2 (2nd column) 


(Type, h, 6, 1) 


lXP,e 

h: chamber height 

6: clearance 

SL: pitch 

S - Straight-through 

0 - 2.75 mm 

0 = 0.25 mm 

0 = 5 mm 

F = Full- interlocking 

1 = 6.25 mm 

2 =0.5 mm 

1 = 8 mm 


3 = 6 mm 


2 = 4 mm 


As figure 8 Illustrates, the calculated values are somewhat lower than the data 
(about 19% for series (s, 0,0,0), 5% only for series (s, 1,0,0)). The trends of the 
calculation are in agreement with those observed in the tests. In particular, K 
is seen in figure 9 to be approximately proportional to inlet swirl and to overaO 
pressure ratio, both in the tests and in the calculations (although, as mentioned, 
with a somewhat lower proportionality factor in the latter case) . 

4.4 Stiffness Coefficients with Shaft Rotation 


Results similar to those in the previous section, but including various shaft 
speeds are given in Table 3 and figure 10. The eccentricity is still static (no 
whirl, Q = 0), and is 0.25 mm in all cases. The parameter E* was used in Ref. 11 to 
correlate entry swirl, and is 
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(17) 


The comparison of data and theory shown in figure 10 indicates more scatter, but 
less systematic deviation than in the cases without shaft rotation (figures 8 and 9) . 
The agreement is best for all the cases with 17 chambers (solid S 3 anbols in figure 10), 
which show an average error of 8.5% and little scatter. 

4.5 Discussion 


The two principal sources of uncertainty in our calculations are the friction 
factors (X*, A”) and the carry-over coefficient B. The friction factor could in 
principle be substantially increased by the relative rotation of shaft and easing, 
since the fluid in each chamber is strongly sheared and develops marked secondary 
flow patterns, leading to enhanced mixing. Examination of data^%or the somewhat re- 
lated case of turbulent pipe flow with swirl does indeed show friction increases of 
up to a factor of four at high swirl. An accurate prediction of wall friction under 
the complex flow conditions of a labyrinth gland is not possible at this time, and 
this is an area requiring more experimental and analytical work. 

The impact of friction factor inaccuracies on calculated cross-spring coeffici- 
ents could be important, although not easily generalizable. In general, the cross- 
forces increase with the deviation between the swirl velocity and its frictional 
equilibrium value. An increase in wall friction in a non-rotating seal will acceler- 
ate the approach to this ultimate swirl, thus reducing the number of chambers where 
the excess or defect swirl is strong and thereby reducing the magnitude of Kyx 
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(whether positive or negative) . This trend is always apparent at large enough fric-* 
tion coefficients, and especially for long multichamber seals. However, for short 
seals and weak frictional coupling a different effect dominates, namely the cross- 
stiffness tends to a limiting value independent of shaft rotation and in the direc- 
tion of the inlet swirl. An example of this behavior for a one-chamber seal is pre- 
sented in Appendix 3, figure A2, where it can be seen that for this particular case, 
increasing friction would lead to increases in |k |. The behavior typical of long 
seals and large friction is illustrated in figure^^ 11, corresponding to case 1 of 
Table 2; here, an increase of the friction factor above the nominal value leads to a 
reduced K , although a reduction by more than about 0.6 would lead to the same 


effect. 


yx 


The carry-over coefficient 3 is clearly another source of uncertainty in the 
model. The sensitivity of calculated cross-stiffness to 3 (or p = 3c ) is illustra- 
ted in figure 12, corresponding to parametric variations on Run 11 of Table 2, and 
figure 13 (from Run 7 of Table 3) . The opposite trend in these two cases is due to 
the fact that in figure 12 the entry swirl is greater than its asymptotic value 
(reached after an infinite number of chambers), while the opposite is true in figure 
13. In both cases, an increase of ]i increases the flow rate q, which has the effect 
of delaying the transition towards the asymptotic c* ; in figure 12 this means higher 
c* in the first 10-12 chambers, with correspondingly larger cross-forces; in figure 
li , the same delaying effect at higher ]X implies lower c* values in the first 6-8 
chambers, and, consequently, lower cross-forces. 

Unfortunately, the state of the art in a priori predictions of 3 is not satis- 
factory. A fuller discussion of this point is given in reference 15. Basically, the 
best-known models (Vermes (ref. 17), Egli (ref. 20), Komotori (ref. 21)), indicate 
substantially different variations of the carry-over factor with number of chambers, 
with - Vermes’ model taking no account of this number at all. The best hope here lies 
with the numerical methods which are now beginning to be applied to internal flow 
problems in seals, although the somewhat primitive state of affairs with respect to 
calculations of fully separated turbulent flows still indicates a need for improve- 
ments. Thus Wadia and Booth (ref. 22) analyzed seal flows with no rotation and ob- 
served discrepancies of up to 13% in calculated flow coefficients when compared to 
data. For dynamics studies in seals, these 2-D or 3-D methods may, in any case, be 
too laborious; their proper role should probably be in furnishing improved semi- 
empirical results for integration into a simple multi-chamber lumped-parameter model, 
of the type considered here. 

5. Parametric Studies 


Reference 15 includes a variety of calculations that illustrate the trends of 
the force coefficients versus variations of seal parameters. Only some of the sal- 
ient results will be mentioned here. 

(a) K increases linearly (but not proportionally) with entry swirl velocity. For 
conditions where the entry swirl exceeds the asymptotic azimuthal velocity, K is 
generally positive, leading to excitation of forward whirl (with respect to the 
swirl direction) . The reverse may be true at lower entry swirls . There is in some 
cases a value of entry swirl at which is zero. 

(b) For multichamber seals with low entry swirl, the first few chambers contribute 

negative K values , while those towards the end of the seal contribute positive 
values. Thus as the seal Is made longer, the sign of may at some point 
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reverse. For seals where the entry swirl exceeds the ultimate azimuthal velocity, no 
such reversal occurs. 

(c) The damping coefficient C , which, together with K controls the side force F , 
is positive in all cases studied, leading to stabilizing^forces of the same order as^ 
those due to This point will be more fully discussed in what follows. 

6. Considerations on Fluid Damping 

An example of calculated damping coefficients is shown in figure 14 (correspond- 
ing to the seal configuration of Run 1 in Table 3). A whirl (critical) frequency of 
739 rad/sec was assumed; at the commonly found ratio 0)/Q = 2 for instability onset, 

= 220 N sec/m , giving = 1.626xlO^N/m. This is several times larger than 

K , and indicates that seal damping forces are in this case sufficient to en^ 

sure stability. Another example of this behavior is shown by the single-cavity 

seal of Appendix 3; here -I- QC is roughly comparable to K (but smaller). 

XX yx 

A very simplified dynamic model will help to put in perspective the roles of the 
different coefficients in stability analysis. Assuming a shaft with mass M and struc- 
tural stiffness K , the equations of motion for small side displacements x,y can be 
combined into 


A . 

M z + C z 4* K z = 0 


(18) 


where C = C + i C , K = K + K + i K and z - x + iy. 

XX xy o XX xy . 

Assuming K , K , /4M and are all small compared to K , as is likely to be 
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the case in practice, we can define the (small) nondimens ional parameters 
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(19) 


the shaft complex 


( 20 ) 


Thus any nonzero k will be destabilizing (in one or the other whirl direction), 
while a negative ^ wifi be always destabilizing; k and g will simply modify 
the shaft natural ^^frequency. Also, the effects of equal vaf2es of |k | and 
- Q C are seen to be equivalent. 

XX ^ 


This discussion has served to indicate that knowledge of the damping factor C 

is at least as essential to studies of fluid-induced destabilizing forces as is 

the side force factor K . Yet, due to the more difficult experimental conditions, 

much fewer data are avaflable on G than on K . 

XX xy 
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We are now in the design stage of a test rig intended to address this problem. 

The general size and flow parameters will be similar to those used by Benckert and 
Wachter, but the sealed shaft will be made to execute forced whirling motion at 
speeds controlled separately from the spinning motion. Pressure distributions will 
be dynamically measured and integrated to produce yalues of the direct and transverse 
forces Fj , F for a range of whirl speeds Q. Both K and C coefficients can then be 
extracted by Ihe methods described in this paper. These features are similar to those 
of a rig described in reference 23 for tests in water. 

7. Conclusions 


A linear analytical model for the prediction of fluid forces in labyrinth seals 
has been presented and discussed. Coiaparison to literature test data shows reason- 
able agreement for the important cross-stiffness K__. The Importance of the damping 


factor has been highlighted and the need for 


yx 


damping data made clear. 
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Squaring Eq. 


Appendix 1. The Unperturbed Solution 
(1) and adding for all chambers yields for the nominal flow rate 
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Also, adding for the first n chambers only gives 
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The momentum equation (Eq. (3)) becomes In the steady state 
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(A-5) 


and X (the Darcy friction factor) given by a modified pipe-flow expression 


, 0.3164 r, . ^ „ 0.25,°hs^^^n 

^ “ 71:25 \ ^2R'^ ^ ^ 

R s 


(A-6) 


Here R is the Reynolds number based on chamber height and the corresponding relative 
flow velocity for fixed or rotating surfaces. The factor sign (v___ ) is needed 
to give the forces , xV their proper direction. Thus, we append the factor sign 
(c*) to X^ and the factor^sign (toR. - c*) to A” , both here and in the first order 
calculations of Appendix 2. is^ the gland hydraulic diameter. Eqs. (A-3) 

through (A-6) can be solved for the distribution c. of azimuthal velocities. In 
particular, the asymptotic velocity (c* ) follows from (A-3) when c| = 
assumed. 


Appendix 2. Coefficients for the Perturbation Equations 


Let the transverse area of a gland (fig. 1) be f. = (h. 4* 5^)£^ in the centered 
position. An asterisk on any variable denotes the undisturbed (centered) condition. 
We obtain for Eq. (7) the following coefficients. 


A = — 
i Y 


* * 
c f. 

^i R Y 

s 


C. 

1 


f 


k 

i 


s 


°1 = 


41 ^1 


2 


kr 

1+1 


* * 
p. q 


* 






y 


2 

1 
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F. = 


P ^ u ^ ^ 

* * 




*1- 


* n 
Ci 


K. = 


* 

f . 

X 

Y 


* 


* 

f , 

M. = ^ ( 


R T 

-J_ 


Y 


c . 

X 


N. = 


* 

2e^ fi 


°1 = 


p *2 2 5*2 

1+1 ^1+1 i+1 

* * 

p. q 


p. = 


P* 5*2 

^1 ^ 1+1 ° 1+1 ^1 " 1-1 ^1 


« 75 O (x" 9 

P. c, , p? 6, 2 


* i 

Cl q 


* 

'1 


8 Y 


U" X" 


* (WR - c.)' 

^ S X 


8c . y 

X ^ 


Q. = 


p *2 * ^2 5*2 

1 - 11-1 ^1 1 
* * * 

Pi q 


^ ^ q . X'U' * , X" U" *s 

\ - — + ~ir~ + "~ir~ ■ ^1^ 

Pi 


s, = 


c. . * 

_izi 3^ 
* * 


„ .3_ (-!_ _:ari_L 

i * ^ * * 

‘’i '* 1+1 =1 '*1 


*1 


* . 

G . £. 

Y = ^ ^ 

i R 

s 


Here y is the ratio of specific heats and R the gas constant. 

s 


Appendix 3. ”Local*^ vs. "Global” Coefficients 


Sample calculations were made for a single-chamber straight-through seal with 
6 = 0.25 mm., il - 8 mm., h = 3 mm., R = 15 cm., 0) - 1000 rad/sec, c = 100 m/sec, 

P - 1.5 atm, P = 1 atm . For a range of whirl frequencies from Q = 0 to Q = 750 
rad/sec, the resulting direct and quadrature forces are shown in figure Al. The 
quadrature force F , which is the one of importance for stability considerations, is 
seen to be very ^ nearly linear with fi, indicating no difference between the local 
and the global definitions of the coefficients (Eqs. (15) and (16), respectively). 
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There is, on the other hand, a slight curvature in the line. The yalues calcu- 
lated for = 500 rad/sec are as follows: 




K Wm) 

K (N/m) 

C (N sec/m) 

C (N sec/m) 



XX 

yx 

XX 

yx 

Local 

(Eqs. (15)) 

610 

11350 

+ 13.7 

- 0.95 

Global 

(Eqs. (16)) 

652 

11480 

13.98 

- 1.02 


For the same seal, with = 0 throughout, figure A2 shows the effect of para- 
metric variations of the friction coefficients (A^ and A” varied simultaneously) at 
various shaft rotation speeds. 

Table 2. K Calculated vs. Experimental (o> ^ ” 0) 

yx t' ^ ^ 

Zx 

C. Exper. Calcula. 

P /p o 5 

Run // Seal Type K a.' o (m/s) U (x 10 N/m) 


1 

S,0,0,0 

17 

.66 

38.4 

.92 

.75 

.611 

2 

S, 0,0,0 

17 

.32 

52.9 

.92 

2.57 

2.091 

3 

S, 1,0,0 

17 

.32 

68.3 

1.02 

1.57 

1.45 

4 

S, 1,0,0 

17 

.66 

33.4 

1.02 

.27 

0.274 

5 

S, 1,0,0 

17 

.66 

48 .5 

1.02 

.423 

0.411 

6 

S,1,0,0 

17 

.56 

39.0 

1.02 

.457 

0.410 

7 

S, 1,0,0 

17 

.79 

38.2 

1.02 

.218 

0.224 

8 

s, 0 , 0,0 

17 

.49 

63.5 

.92 

1.89 

1.586 

9 

S,0,0,0 

17 

.39 

54.3 

.92 

2.22 

1.747 

10 

S, 0,0,0 

17 

.79 

15.6 

.92 

.184 

0.160 

11 

S, 0,0,0 

17 

.49 

64.2 

.92 

1.75 

1.605 

12 

S, 0,0,0 

17 

.49- 

34.5 

.92 

1.05 

0.807 

13 

s, 1 , 0,0 

17 

.49 

82 . 68 

1.02 

.98 

1.098 

14 

s, 1 , 0,0 

17 

.49 

40.78 

1.02 

.57 

0.509 

15 

S,0,0,0 

17 

.32 

38.2 

.92 

1.9 

1.431 

16 

S, 0,0,0 

17 

.32 

27.6 

.92 

1.2 

0.937 

17 

F.3,2,1 

9 

.49 

144,7 

.665 

1.47 

1.606 

18 

S, 0,0,0 

17 

.32 

45.5 

.92 

2.39 

1.757 

Average error 

(in absolute value) ~ 18 

.3% 



(4. 5% 

for (S,l 

,0,0) , 

18 .6% 

for (S, 0,0,0)). 




Table 3. K Calculated vs. Experimental (w 0) K . 

yx c (ii y- — — — 

yp *Q rot Exper. ^Calcula. 

Run y/ Seal Type K a^ o (m/s) (rad/s) u (x 10 N/m) 


1 

F,3,2,2 

17 

.66 

43.2 

1000 

.66 

.189 

0.177 

2 

F.3,2,2 

23 

.66 

47.1 

1000 

.66 

.44 

0.349 

3 

F.3,2,2 

23 

.66 

40. 

1000 

.66 

.38 

0.315 

4 

F.3,2,2 

23 

.79 


993.3 

.66 

.307 


5 

F,3,2,2 

23 

.793 

66.2 

993.3 

.66 

.323 

0.398 

6 

F,3,2.2 

23 

.793 

49.8 

746.67 

.66 

.189 

0.228 

7 

F.3,2,2 

17 

.793 

50.7 

993.3 

.66 

.248 

0.230 

8 

F,3,2,2 

23 

.657 

86.74 

993.3 

.553^ 

.442 

0.398 

9 

F.3,2,2 

23 

.657 

83.64 

746.67 

.606^ 

.290 

0.340 

10 

F,3,2,2 

17 

.66 

27.9 

1000 

.66 

.112 

0.088 

11^ 

F.3,2,2 

17 

.66 

43.2 

1000 

.66 

.20 

0.187 

12^ 

F.3,2,2 

17 

.66 

27.9 

746.67 

.66 

.12 

0.0849 

13^ 

f.3,2,2 

17 

.66 

15.5 

500. 

.66 

.044 

0.0178 

14" 

f.3,2,2 

23 

.66 

47.1 

1000 

,66 

.40 

0.328 

15" 

F, 3.2.2 

23 

.66 

30. 

746.67 

,66 

.23 

0.170 

16" 

F, 3,2.2 

23 

.66 

54.7 

500 

.66 

.09 

0.123 


j. — jji calculated from measured mass flow rate: friction factor for a 
channel used. 

2 -- resultant forces from chambers 7 - 17 only, 

3 — resultant forces from chambers 7-23 only. 

Average error (in absolute value) “ 23.0% 
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Figure 6. “ Calculated versus experimental axial pressure distribution 
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Figure 7. - Calculated versus experimental circumferential velocity 
distribution. 
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Figure 8. - Calculated versus experimental cross spring coefficient Kyx 
(“rot “ • 
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Figure 9. - Entry swirl parameter and calculated and experimental cross spring 
coefficient Kyx. 














Figure Al. - Direct and quadrature forces in a short seal. 



Figure A2. - Effect of friction factor variation on cross-stif fness in a short 
seal. Inlet swirl velocity, Gq, 100 m/sec. 
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